Chronotopic Lyapunov Analysis: 
(I) a Comprehensive Characterization of ID Systems 



Stefano Lepri 

Dipartimento di Fisica, Universitd di Bologna and Istituto Nazionale di Fisica Nucleare 

1-40127 Bologna, Italy 

Antonio Politi 

Istituto Nazionale di Ottica and Istituto Nazionale di Fisica Nucleare 

1-50125 Firenze, Italy 

Alessandro Torcini 
Theoretische Physik, Bergische Universitdt-Gesamthochshule Wuppertal 
D-42097 Wuppertal, Germany 
(February 5, 2008) 



Abstract 



Instabilities in ID spatially extended systems are studied with the aid of 
both temporal and spatial Lyapunov exponents. A suitable representation of 
the spectra allows a compact description of all the possible disturbances in 
tangent space. The analysis is carried out for chaotic and periodic spatiotem- 
poral patterns. Singularities of the spectra and localization properties of the 
associated Lyapunov vectors are discussed. 
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I. INTRODUCTION 



Linear stability analysis of chaotic dynamics is usually concerned with the problem of 
measuring the divergence of nearby trajectories [|l|. In extended systems, the spatial de- 
pendence of the state variable requires considering also propagation phenomena and, more 
generally, spatial inhomogeneities 0. The full characterization of a generic perturbation 
involves both its temporal and spatial growth rates as complementary measures of its in- 
stability properties. Two classes of Lyapunov exponents have been separately introduced 
for this purpose: the former aims at describing the temporal evolution of disturbances with 
an exponential profile in space 0; the latter deals with the spatial shape of a perturbation 
defined on a given site at all times 0]. These, which will be called as temporal and spa- 
tial exponents, represent the starting point of this paper. Their definitions are reported in 
section III. 

Other indicators have been introduced in the attempt of characterizing the mechanisms 
of information transport. In this case, the object of study is the propagation of initially 
localized disturbances and it led to the introduction of the comoving Lyapunov exponents 
I^J^. Moreover, the Lyapunov analysis of a given (H-l)D pattern can, in principle, be 
carried out along any direction. In fact, by performing a spatio-temporal rotation, i.e. by 
combining the role of space and time, it is possible to define still another class of "rotated" 
exponents 0. 

It has been already shown that the maximum comoving Lyapunov exponent can be ob- 
tained from the maximum temporal exponent through a Legendre transform 0. Moreover, 
the Kolmogorov-Sinai entropy turns out to be independent of the propagation direction 
along a spatio-temporal pattern 0. Therefore, there must exist strong relationships among 
the several classes of exponents so far introduced. A particularly appeahng problem is that 
of identifying the independent indicators which are necessary and sufficient for a complete 
characterization of a linear instability. 

The present work, divided into two parts, intends to represent a first step towards a 
general and coherent theory. In this first part, we start from the observation that the 
most general perturbation is identified by two rates A, /i, describing its growth in time and 
space, respectively. The number of linearly independent perturbations can be expressed 
by two integrated densities, nx{fi, X) and A). They provide a dual representation of 

the problem since any perturbation is unambiguously identified by the pair {n^,n\) as well 
as by (/i. A). The properties of these indicators are studied in several models of coupled 
map lattices (CML), since we are confident that the tools and the results apply equally well 
to systems with continuous space and time variables. The relationship with the indicators 
arising from reference-frame-dependent analysis of disturbances will be thoroughly discussed 
in the second part p|. 

In section II, several classes of CML-models are recalled, which have been designed 
to mimick reaction-diffusion systems, many-degrees-of-freedom Hamiltonian dynamics, and 
systems with a conserved order parameter. As already anticipated, in section III spatial 
and temporal spectra are formally defined. In section IV, we introduce the so-called (/x, A)- 
plane, which provides a compact graphical representation of Lyapunov spectra. It also allows 
enlightening analogies and differences between spatial and temporal exponents. 

CML Lyapunov analysis has several analogies with the Schrodinger problem in disordered 
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lattice systems [1,0. This allows extracting from the spatial spectra information about the 
localization of the temporal Lyapunov vectors, thus confirming the close relationship between 
the two approaches. A discussion of the analogies with the Anderson tight-binding problems 
is presented in section V, where the the {fi, A) representation is used for a straightforward 
identification of exponentially localized vectors. Moreover, the presence of of singularities 
in the spatial Lyapunov spectrum is interpreted as a stretched-exponential behaviour of 
the vectors. Finally, another type of singularity occurring in the temporal spectrum is 
illustrated. Some general remarks introducing the reader to the second Part are presented 
in Section VI. 



II. ID MODELS FOR EXTENDED SYSTEMS 

Spatio-temporal chaos and instabilities in extended systems have been widely studied 
with the aid of simplified models of react ion- diffusion processes, whose ID form is of the 

type i 



dty = Fiy) + Ddly 



(1) 



with the state variable y{x,t) defined on the domain [0, L] (periodic boundary conditions 
y{0,t) = y{L,t) are generally assumed). The nonlinear function F accounts for the local 
reaction dynamics, while the diffusion constant D represents the strength of the spatial 
coupling. 

Unfortunately, accurate numerical investigations of partial differential equations like (|Tp 
can require very long CPU times even on powerful computers. The introduction of simplified 
models has been of great help for understanding the statistical properties of spatio-temporal 
chaos. A first step in simplifying model (|lD is achieved by discretizing the space variable, 
i.e. by considering coupled oscillators on a lattice. This procedure is generally justified 
whenever the fluctuations are cut off below a certain spatial scale. After replacing the time 
derivative with its finite difference approximation, one is led to the CML dynamics of the 
form 



yU, = F{y\) + (1 



Vn 



(2) 



with n being respectively space and time indexes labeling each variable y]^ of a lattice of 
length L and e = 2D gauges the diffusion strength. Model (^) has some technical limitations 
since, for arbitrary F, it leads to diverging solutions. To avoid such difficulties one can 
postulate a CML dynamics of the form |[To|JTT|| 



y 



n+l 



Vr. 



(3) 



where / is a nonlinear function mapping a given interval / of the real axis onto itself. 
Periodic boundary conditions yl^^ = y]^ are again assumed. A generalization of model @ 
has been proposed p to study convective instabilities, namely 



y^+i = /((!- e)yl + s[{l- a)yl' + ayl^']) 



(4) 
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The parameter a (comprised between and 1) accounts for the possibihty of an asymmet- 
ric couphng, mimicking open-flow systems (first order derivatives in the continuum hmit). 
Obviously, the symmetric case (^ is recovered for a = 1/2. 

Coupled maps are sometimes used also as test ground for extended Hamiltonian systems. 
The corresponding models are generally formulated in such a way that the evolution in 
tangen space is described by symplectic matrices. As a specific example we will refer to the 
coupled-standard-maps lattice 



In + Pn+1 



sin 9n + I sin(g^^^ - O + | sin(g; 



(mod 27r) 



(5) 



where the parameter k rules the amplitude of the nonlinear coupling. 

Finally, there is a last class of models that has been used to study evolution in the 
presence of a conserved order parameter [|1^]. Such CML models have been introduced 
in order to make a closer contact with the dynamics of microscopic conserved dynamical 
variables in strongly turbulent fluids. A general scheme for constructing a CML with a 
conserved quantity is to impose that Hn^i — Un is written as a spatial gradient. For example 



Vn+i = yn + 9{ylt^) + giiin ') - '^givn) 



(6) 



preserves the quantity X^i Un- 



III. LYAPUNOV ANALYSIS OF SPATIO-TEMPORAL CHAOS 

Lyapunov exponents are the main statistical tool for the study of low- dimensional strange 
attractors, since they measure the linear instability of trajectories in phase space. Moreover 
they can be related to other ergodic indicators such as entropies and dimensions |^. In 
extended systems, whenever many interacting modes are present, both temporal and spa- 
tial instabilities must be taken into account. Two complementary approaches have been 
developed which are based, respectively, on the temporal (spatial) growth rate of space 
(time) periodic perturbations. The former approach is the standard method to determine 
the spectrum of Lyapunov exponents. The latter one, which leads to introduction of spatial 
exponents, is useful both in revealing the localization properties of Lyapunov vectors |^ and 
characterizing the stability of open-flow systems [0. 

Numerical simulations [|l5l and approximate analytical results ||l6l indicate that a well 



defined thermodynamic limit for the set of Lyapunov exponents is reached for large enough 
system size, both in space and time. As a consequence Kaplan- Yorke dimension and 
Kolmogorov-Sinai entropy are extensive quantities so that the dynamics of a long chain 
can be roughly seen as that of many independent adjacent subchains. This fact permits to 
define Lyapunov spectra independent of the system size (see below). 

The following two subsections are dedicated to a brief review of the above mentioned 
definitions. For simplicity we limit ourself to the symmetric CML case Eq. (^), the extension 
to spatio-temporal flows being straightforward. 
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A. Temporal Lyapunov Exponents 



The standard Lyapunov spectrum is obtained by following the evolution of a perturbation 
6yl^ in the tangent space with periodic boundary conditions assumed for it. A more general 
class of Lyapunov exponents has been introduced in Ref. @] by enlarging the set of allowed 
perturbations to include exponentially shaped profiles such as 6yl^ = ^^e'^*. The evolution 
equation for the scaled variable $^ reads as 



(7) 



where = /' (^^Vn ^ + (1 — e)yl^ + fy^"*"^) is the local multiplier, and $^ obeys periodic 
boundary conditions. The exponents Aj(/i) (1 < j < L) obtained from the iteration of Eq. 
(0) for a fixed value of fi will be hereafter called temporal Lyapunov exponents |T^. The 
resulting temporal Lyapunov spectrum (TLS) is defined in the limit L — > oo, as 

A(/i,nA) = XjiiJ,) , nx = j/L (8) 

where the integrated density (the fraction of exponents smaller than A) nx ranges, by con- 
struction, between and 1. For fi = the standard TLS is recovered. Exactly the same 
arguments apply to continuous-time models as, e.g., chains of coupled oscillators. When the 
space too is continuous, the only difference is that nx is unbounded from above. 

It it well known that the computation of Lyapunov exponents from Eq. proceeds 
through the construction of an ordered sequence of Lyapunov vectors (i.e. an orthonor- 
mal basis), which generate the most expanding subspaces of increasing dimension In 
general, the direction of such vectors depends on the point in phase-space where they are 
computed. These fluctuations make an analytic treatment of the linear stability problem 
hardly feasible. However, their spatial structure can be clarified by following in space the 
evolution of perturbations. This is precisely the subject of the next section. 



B. Spatial Lyapunov Exponents 



Spatial Lyapunov exponents have been introduced as a tool for investigating the structure 
of temporal Lyapunov vectors |Q. By assuming a temporal exponential profile for the 
perturbation, 6yl^ = e'^"\E'^, the evolution equation in tangent space reads as 



£^i-i + fl - eW + -^-^^ 



(9) 



which, in turn, can be seen as a spatial recursive equation. 



i+1 



(10) 



with boundary conditions = ©l, ^t+i — ^i- (0) starting point 

of a transfer matrix approach. In fact, by introducing the column vector v(z) = 
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{Q\, Q\, . . . , 9^, "^l, . . . , the spatial evolution can be described in terms of products of 
matrices of the form 

where and 1 are the (T x T) null and identity matrix, and A* is a bidiagonal matrix, 
with = —2(1 — e)/e and J^n.n~v\ = 2e'^/(£:m^). Since det J*=l the spatial dynamics is 
"conservative" . Numerical implementation of this method requires the knowledge of the 
local multipliers m^, which have to be computed, in principle, for a generic orbit of period 
T. In practice, it is assumed that boundary conditions are irrelevant in the thermodynamic 
limit T — > cxD, and that any initial condition chosen according to the invariant measure is 
almost equivalent. 

For a given T, there exist 2T Lyapunov exponents /ij. This doubling in the number of 
the degrees of freedom is related to the two-step spatial memory of Eqs.(|IO|). Moreover, 
being J* = the spatial Lyapunov spectrum (SLS) associated with any symmetric 

trajectory is invariant under "space" -reversal i —i, i.e. fij = —fiT+i-j- Notice that this is 
no longer true for models like (H) with a 7^ 1/2. In analogy to the TLS, the SLS is defined 
in the limit T ^ 00, as 

j — 1/2 

fi{X,n^) =fij{X) , = — 1 , (12) 

where j denotes the j-th spatial Lyapunov exponent. The term 1/2T — 1 is a finite-size 
correction, introduced to fix the center of symmetry at = 0, independently of the size T. 
The density of spatial exponents ranges between -1 and 1. 

In analogy with disordered systems, the minimum positive exponent fi{X, 0) can be in- 
terpreted as the inverse of the localization length / of the Lyapunov vector, provided that A 
belongs to the TLS. 

The direct implementation of the above method faces the difficulty that iteration of the 
nonlinear equation in space usually does not converge onto the same attractor as obtained by 
iterating in time the original model. In fact, it turns out that the invariant spatio-temporal 
measure corresponds to a strange repeller in space (this point will be elucidated in the second 
part). Thus, one must first generate a 2D pattern of local multipliers and then use it in the 
computation of the spatial Lyapunov exponents. 



IV. THE (/i, A)-SPACE 

A more symmetric representation of both spatial and temporal spectra is obtained by 
expressing Eqs. (H) and (|I^) as n\{fi, A) and n^(/i. A). Both functions are well defined in so 
far as A and are monotonously decreasing for increasing density. Any point (/x. A) identifies 
a specific perturbation growing as e^" in time and as e^* in space. The integrated density 
(nx) corresponds to the normalized number of "modes" with spatial (temporal) exponent 
smaller than fi {X). 

The set V of points in the (/i, A)-plane corresponding to admissible perturbations is 
identified by simultaneously requiring. 
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(13) 



One should notice a similarity with the dispersion relation for wave propagation: given a 
specific spatial profile, the temporal growth rate is bounded in a given interval from the 
above inequalities. The boundaries and dT>x of the domain T> are obtained by setting 
the l.h.s. equal to zero in Eq. ([TB|). In Fig. 1 we report such borders for a chain of 
coupled logistic maps (/(x) = 4a:(l — x)) and e = 1/3 (a), e = 2/3 (b). As the two borders 
coincide in both cases we have a first indication that and nx are not independent of one 
another. Although this is not surprising, since both quantities arise from the iteration of 
the same linear set of equations, it is far from obvious. The invariance of the dynamical 
equations under the transformation i —>■ —i implies a left-right symmetry of the border dV 
(we drop the subscript whenever dV^j, = dV\). The upper and lower boundaries at /i = 
correspond to the standard maximum Xmax and minimum Xmin Lyapunov exponent (for 
e = 2/3, Xmin = — oo 0). Three typical representatives of temporal and spatial spectra are 
reported in Fig. 2 (a-c) and (d-f), respectively, with reference to e = 1/3. The SLS obtained 
for Xmin < A < Xmax IS characterized by a single band, while for larger and smaller values 
of A two symmetric bands appear. The TLS appears to be characterized by a single band 
for any value of /x. For e = 1/3 the size of the band increases monotonously up to /i = fic, 
where the minimum exponent diverges to — oo, while in the latter case the divergence occurs 
at /i = 0. Above fic, the size of the band shrinks to zero and its position increases linearly 
with /i. 

This can be easily understood by realizing that, for large /x, Eq. (|^ reads as 

K+i^^n^e^K^' . (14) 

Accordingly, the asymptotic value of the Lyapunov exponent is 

A^/i + log(£/2) + (log|m^|) , (15) 



where the average (■) should be taken along the line i = I — n for any given /. Eq. (|T5D 
is increasingly accurate for fi —>■ oo, when the evolutions along the lines i = I — n are 
exactly decoupled, so that an infinite degeneracy in the spectrum is necessarily found. It is 
transparent from Eq. (^) that the slope of the branch must be equal to 1. Such a value is 
nothing but the limit velocity for the propagation of disturbances in a lattice with nearest- 
neighbour interactions. In fact, in the case of a coupling extended to s neighbours, Eq. 
(|1^) becomes, $^4.1 ~ m^ee^^$5^"^^/2, so that the slope is, in general, equal to the number 



of neighbours s. Finally, notice that log(e/2) + (log|m^|) indicates the intersection of the 
asymptote with the vertical axis. 

The above described phenomenology turns out to be rather general: simulations per- 
formed with several different maps exhibit the same behaviour as long as one limits the 
analysis to nearest neighbour coupling. 

An equivalent but instructive representation is obtained considering nx and as inde- 
pendent variables. In fact, it turns out that any pair {n^,nx) identifies unambiguously a 
spatio-temporal perturbation, exactly in the same way as /i and A do. The two integrated 
densities are, by construction, normalized in the rectangle [— 1,+1] x [0,+l]. Numerical 
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simulations performed on logistic maps show that the domain V in the (/i, A)-plane corre- 
sponds to a domain B which, strange enough, does not coincide with the full rectangle and 
indeed only points with nx > are obtained (see Fig. 3). We have no explanation for 
this phenomenon. 

The four lines identifying the border dT> are mapped onto the four distinct points P = 
(0, 0), Q = (0, 1), R+ = (1, 1) and R- = (—1, 1). A detailed characterization of the mapping 
in presented in Fig. 3. Symmetry reasons require that the line = is mapped onto the 
line = 0. All vertical lines in the {fi, A) plane correspond to curves departing from P in 
{n^,nx). Their ending point E depends on fi: for < —fic, E = _R+, for —^c < < /^c, 
E = Q, and for fi > fi^ E = R_. The critical curves for fi = are straight lines arriving 
at (±0.5, +1). 

The images of all the horizontal lines connect R- with i?+ exhibiting an obvious sym- 
metry with respect to n^. All curves corresponding to A < Xmin pass through the point Q. 
Analogously, those ones with A > Xmax pass through P. 

We can thus conclude by remarking that there is a sort of duality between the two 
representations. In fact, the four segments delimiting B correspond to 4 points in the {fi, A) 
plane, more precisely to (±cxo, +oo) and (±/Xc, — oo) in Fig. 1. 

In order to examine the generality of the above scenario in the following subsections we 
discuss two specific cases: (a) spatially homogeneous and stationary pattern, which can be 
analytically treated, and (b) spatio-temporal periodic orbits. 



A. Homogeneous Chains 

The simplest chaotic dynamics one can think of is the evolution of piecewise linear maps 
of the type f{x) = rx (mod 1). In this case, the local multiplier is everywhere constant 
(both in space and time). The resulting Lyapunov spectra coincide with those of spatio- 
temporal fixed points (stationary and homogeneous solutions). The TLS is easily obtained 
by noticing that the associated eigenvectors are nothing but the Fourier modes e*''-' (where 
j is the imaginary unit) of the chain. The TLS can be determined by Fourier transforming 
Eq. d^), which allows calculating the multiplier m for a given wavenumber q. The Lyapunov 
exponent X{fi,n\) = log|m| is 



X{n,nx) = logr + - log 



^1 — e)^ + 2£(1 —£:) cosh /i cos TTTT-A + £^(cOsh^/i — sin^TTTT-A) , (16) 



where nx is equal to the spectral density g/vr, since in virtue of the node theorem, all the 
multipliers are naturally ordered with q. 

The same technique can be used to evaluate the SLS, the only difference being the two- 
step spatial memory in Eq. (0). The spectral problem can be solved by determining the 
eigenvalues of suitable 2x2 matrices, obtaining 



fj.{X,n^ 



1 



cosh 



-1 



;2)2 + e'^e'^^ sin^fvm, 



(17) 



where 
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r^{l — e) — 2r(l — e) cosTin^e 



2A 



(18) 



Eqs. (|T6|,|T7D can be inverted to give the integrated densities 



± 1 
= — arccos 

^ TT 



— r(l — e) coshyU ± Yr2(l — 2£:) sinh fi + e 
er 



2A 



(19) 



n„ = — arccos 



-r(l — e) sinh ± cosh/iA/r2(l — 2£:) sinh fi + e 



o2\ 



(20) 



For e > 1/2 



For e < 1/2, only the positive solution exists so that nx = and = 
both solutions are acceptable. In fact, the spectrum of multipliers m passes through zero 
and then n(log|m|) is no longer invertible. In this case , two branches exist, ordered in an 
opposite way with respect to A. The integrated density of Lyapunov exponents is obtained 
by summing the contributions arising from the two partly overlapping bands. 



and 



n„ 



n 



n 



n. 



n, 



for -1 < < 
for < < 1 



The analytic expressions of the border &D can be explicitely obtained by imposing that the 
argument of the arccosinus function in n\ is ±1. For e < 1/2 this procedure yields 



A = logr + log \e cosh /x ± (1 — e) 



(21) 



For e > 1/2, the upper border is unchanged (plus sign), while the lower one is now obtained 
from the identification of the inversion point where nj 



n 



A ) 



A = log r + log 



{26 — 1) sinh n 



(22) 



It should be noticed that the expressions of the borders, as obtained from Eq. ( pj]) or Eq. 
(^) are identical. The two borders reported in Fig. 1 for logistic maps are qualitatively the 
same as those obtained in the present case for e < 1/2 and e > 1/2, respectively. 
Moreover, from Eqs. (^,^), 



±cosh~^(i^) ife<l/2 
otherwise 



The value of /ic is independent of the map /. In fact, the Jacobian matrix can be factorized 
as the product of two matrices, the former containing the detail of the local dynamics, the 
latter being the diffusion operator, which depends only on /i and e. It turns out that it is 
precisely the determinant of the latter operator which vanishes at /i = fid^)- This singular 
behaviour is associated with the existence of a Fourier mode, (—1)*, which is killed in just 
one iterate. 
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Extension of the above analysis to asymmetric CML (^) yields the following expressions 
for the borders 



A =logr + log|ee + 2a sinh/i ± (1 — e)| 

(1 - a)e-'^) 



A = log r + lo 



4a(l— a) 



(23) 
(24) 



which generalize Eqs. (|2l| , |22|) , respectively. Two examples of boundaries are reported in 
Fig. 4 for a = 3/4 and a = 1. The curves are still symmetric, though now the center of 
symmetry is /t = log[(l — a)/a]/2. For decreasing (increasing) values of a from 1/2, the left 
singularity moves towards — oo (+C)o) and eventually disappears at a = (1). 

Finally, one should notice that all the results contained in this subsection apply also 
to the case of a homogenous random evolution: the only difference is that logr is to be 
everywhere replaced with the average value (log \ml^\) . 



B. Periodic Orbits 



An interesting representation of the stability properties in the (/i, A)-plane is obtained 
for periodic (both in space and time) orbits. It was already noticed that the standard 
TLS exhibits a band structure. This feature generalizes to the presence of a series of holes 
within the domain T>. The structure of the temporal spectra is determined from the spatial 
periodicity and vice versa. More precisely, the number of bands of the TLS (SLS) can be at 
most equal to the spatial period (twice the temporal period). The borders &D are reported 
in Fig. 5 for a hypothetical orbit of period 3 in space and 2 in time. Since the reference 
to a specific CML model is not of particular relevance in this case, the values of the 6 
multipliers have been chosen a priori without any relation with a preassigned dynamics. As 
the integrated densities remain constant inside each hole, the latter ones are mapped onto 
single points in the {nx,n^) plane which are analogous to P, Q and R±. 



V. LOCALIZATION IN TANGENT SPACE 

The localization of Lyapunov vectors observed in Refs. P,pO|, implies that a generic per- 
turbation is sensibly amplified in a limited spatial region. Heuristically this phenomenon can 
be explained by noticing the analogy with the Anderson localization problem. The presence 
of the discrete Laplacian operator makes Eq. (^ similar to a wave equation on a lattice, or 
to a Schrodinger equation in the tight-binding approximation. Chaotic fluctuations of the 
multipliers play the same role as a time-dependent disorder in the above problems. 

A further analogy exists between the spatial evolution and the stationary 2D Anderson 
model, defined by the equation 

+ + ^'n-l + C+1 = - V:Wn (25) 

where cu is the rescaled energy of the eigenstate V'n is the random potential bounded 

between —W/2 and 1^/2. In this case, both u and W can be interpreted as control param- 
eters in a dynamical equation. The main difference with Eq. (|T0]) is that now also the index 
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n denotes a spatial direction. The lack of time-reversal invariance in the CML is reflected 
by the absence of a term proportional to ^n_i- 

On the basis of this simple analogy, one might expect that spatially chaotic solutions of 
Eq. (|10D , characterized by pseudo-random sequences of multipliers m^, are associated with 
localized vectors. In fact, the theoretical analysis of the 2D Anderson model predicts the 
existence of exponentially localized eigenf unctions for arbitrarily weak disorder . On the 
other hand, numerical simulations show that /x(A, 0) is always for A belonging to the 
standard TLS, implying that the localization is not purely exponential. This contradiction 
can be resolved by noticing that in the CML case, at variance with the Schrodinger problem, 
the transfer matrices are iterated assuming an exponential profile along the transverse di- 
rection. Accordingly, different operators are involved in the two procedures, with an evident 
symmetry breaking in the latter case. 

Since yu(A, 0) yields the localization length /, it is reasonable to speculate that, whenever 
/i(A, 0) = 0, the behaviour of the SLS in the vicinity of = gives information about a 
supposedly lower-than-exponential decay of the Lyapunov vector. A first necessary condition 
for this conjecture to be true is that the convergence (for T oo) of the SLS to its 
asymptotic shape is sufficiently fast so as to guarantee that the finite-size value of the 
minimum positive exponent is essentially determined by the difference of the corresponding 

from 0, i.e. 

/^t(A)^/.(a,-^) . (26) 

This point will be discussed in the last section of this section; here we limit ourselves to 
notice that numerical simulations reveal that the convergence is not always sufficiently fast. 
Assuming a power-law behaviour of the SLS, i.e. 

/i ~ |n^.|^ n^^O , (27) 

and Eq. (^) to hold, then the vector 6yl^ decays as exp[— L//(T)] ~ exp[—L/T^']. Whenever 
there is a perfect symmetry between spatial and temporal directions, as in the Anderson 
model (|25|) , T can be replaced by L in the above expression yielding a decay behaviour 

6y^; ~ exp[-L(i-^)] . (28) 

Accordingly, if < 7 < 1, a stretched-exponential envelope is obtained. We should anyhow 
recall that T and L are not in general interchangeable, so that the above picture is only 
approximately correct. 

A. Frozen Random Patterns 

For stationary but spatially chaotic states, the evolution in tangent space is mapped 
exactly onto a ID localization problem. In this case, it is rigorously proved that all the 
eigenvectors are exponentially localized. Thus, qualitative differences with the scenario 
depicted in the previous section are expected. Indeed, frozen random patterns are, to our 
knowledge, the only case where the borders derived from spatial and temporal spectra do 
not match. This is illustrated in Fig. 6, where it is clearly seen that the spatial Lyapunov 
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exponents p2[ are all bounded away from zero. As a consequence there is a pair of borders 
which cannot be revealed by the temporal analysis (see dashed line in Fig. 6). In fact, upon 
increasing fi from 0, both the maximum and minimum temporal Lyapunov exponent are 
constant until the spatial border is reached (for clarity reasons, only the behaviour of the 
maximum exponent is reported in Fig. 6). At larger values of /i, the two borders coincide. 

The anomalous structure of "D is a consequence of the exponential localization of the 
temporal Lyapunov vectors. In fact, the dashed curve is nothing but l/l for all temporal 
exponents belonging to the spectrum. 

A further interesting feature concerns the identification of the critical value fii, where the 
upper temporal border departs from the spatial one. By noticing that the upper border is the 
Legendre transform of the maximal comoving Lyapunov exponent P|, it appears natural 
to conjecture that &Dx can be obtained from dV^^ by means of a standard construction to 
remove all changes of concavity. In other words we expect /x; to be located at the inflexion 
point of &D^. This seems indeed to be true. We shall comment more extensively about this 
point in the second part of the paper. 



B. Symplectic Maps 

There is a strict analogy between the application of the transfer matrix approach to the 
Anderson model ( p5l) and the computation of Lyapunov spectra in Hamiltonian systems, 
since in both cases one deals with products of "random" symplectic matrices [^. Generally 
speaking, both the symplectic structure and the time-reversal invariance induce a "pair- 
ing" in the Lyapunov spectrum which turns out to be invariant under the transformation 
A —A. Consequently, a similar structure of the borders is expected along the and 
A directions. Obviously, in the Anderson model both are spatial directions and they can 
be exactly interchanged, so that we expect the exact invariance under the transformation 
(/i. A) (A,/x). 

The main difference with respect to the previous class of dissipative CML is the slope of 
the lower asymptotes, which are now tilted at a ±7r/4 angle. This is precisely the consequence 
of the invertibility of mapping (|^). Indeed, when the evolution is reversed in time, the 
maximum velocity of propagation of a disturbance is 1 (when nearest-neighbour coupling is 
assumed). This fact allows to interpret the vertical asymptotes in Fig. 1 as an indication 
of an infinite (backward) velocity. The non-invertibility of the dynamics leads to a global 
unavoidable ambiguity of the preimage already in one step. 

Simulations of model (|) show two qualitatively different domain structures. Indeed, 
depending on the nonlinearity k, a hole may exist around (A,/i) = (0,0), indicating that 
the standard TLS is made of two distinct bands. The two possibilities are illustrated in 
Fig. 7, where &D is reported for k = 0.5 and k = 4.0. At variance with the latter case, 
at k = 0.5 a substantial disagreement is found between spatial and temporal borders. This 
artifact seems to be caused by strong finite-size effects in the convergence of the SLS to its 
asymptotic form. This may be connected with the slowly decaying temporal correlations 
(up to time of order 10^) previously revealed by the analysis of diffusion properties and 
of effective Lyapunov exponents 



To close this section, we want to remark that the theory developed in [21| for the 2D- 



Anderson problem, predicting that all states are exponentially localized, would imply the 
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presence of a hole at (A, fi) = (0, 0) for any value of k in the present case. On the other hand, 
direct numerical simulations of model (|25| ) appear to be quite inconclusive due to the rapid 
divergence of the localization length for W ^ [Q. Thus, also in our dynamical model 
it is not possible to conclude whether there are two truly distinct regimes as suggested by 
Fig. 7. 



C. Localization and S pat iotemp oral Chaos 

In this subsection we turn again to the general case of spatio-temporal chaos, reporting 
some results on the power-law singularity in the SLS. In general, for A belonging to the 
standard TLS, the corresponding SLS exhibits a linear behaviour around = 0, that is 
7 = 1 (see Eq. (P7|)). Hence, according to Eq. (P^), the Lyapunov vector decays more slowly 
than a stretched exponential. In fact, direct numerical investigations reveal a power-law 
localization 



The only exception to the above picture is represented by the extrema of the standard 
spectrum, i.e. A = Xmin, ^max, where 7-values substantially smaller than 1 are observed. 
This is not surprising in that these are the points where a gap is closing in the SLS (see for 
example Fig. 1). Although this singularity arises at a band merging point in the (/x. A) plane, 
7 is not trivially linked with the shape of the boundary around /z = 0. In fact, while 7 may 
depend on the model, the shape of dV around = turns out to be always parabolic. This 
has been verified both analytically for homogeneous chains by expanding Eq. (|I^ around 
/i = 0, and numerically for various models (see Fig. 8, where 6\ = A(/i, 0) — Xmin is plotted 
versus /i). 

Extensive numerical simulations have been performed on several classes of CMLs to 
measure 7, namely model (^ with logistic, tent, skewed and Bernoulli maps, model with 
linear coupling (^, and symplectic maps (|^). The mapping functions as well as the measured 
values of 7 are reported in Table 1, indicating that the exponent is either approximately 
equal to 1/3 or to 1/2. The inaccuracy is essentially to be attributed to finite-size effects, as 
the maximum temporal exponent has been always computed with a relative error of 10~^. 
Actually, the above estimates can be improved by assuming that the spectrum is an analytic 
function of n^. More precisely, the power-series expansion of to the second order, 

l^^anj^ + bnl'' , (29) 

can be used to fit the numerical data in the vicinity of yU = 0, in order to determine the 
parameters a and b for a given 7. The best estimate of the exponent 7 is finally determined 
as the one which optimizes the fit. As a result of this application the 7-values are even 
closer to 1/3 and 1/2 confirming that there are two universality classes. The analysis of the 
considered models reveals that the exponent 1/2 is obtained only when the local variable is 
scalar and if the local multipliers are all positive (negative). Hence the exponent 1/3 seems 
to arise from "dephasing" effects due to a local rotation of the perturbation (if the variable is 
a vector) or to the randomness of the sign of the multipliers (if the variable is a scalar). The 
larger generality of the latter value is further confirmed by its observation in very different 
contexts such a band random matrices [EB] or nonequilibrium molecular dynamics 



A further interesting interpretation of the two universality classes appears from the 
structure of the temporal evolution in tangent space (Eq. (^). Pikovsky and Kurths have 
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shown that the equation resulting from the change of variable /i^ = In 5?/^ is approximately 



equal to a discrete version of the Kardar-Parisi- Zhang equation |29] for the scalar field h{x, t) 
{x and t replace i, n, respectively) 

dth = r],dlh + V2{d,hf + C, (30) 

where ^{x,t) = Inm^ is a noise term arising from space-time chaos, and 771^2 are two con- 
stants. A first limitation to an exact correspondence stems from higher order spatial deriva- 
tives the presence of which should not affect the scaling behaviour ||2^. A more relevant 



difference arises when the sign of is a fluctuating quantity, since the definition itself of /i^ 
is meaningless, as already argued in According to our previous considerations, 7 turns 
out to be 1/2 if and only if the reduction to Eq. ( pOD is possible. In such cases, from Eq. (p8[), 
we expect for the Lyapunov vector a stretched exponential profile Sy^ ~ exp(— L^/^), which 
is perfectly consistent with the prediction that the field h describes a Brownian motion in 
space. Thus, the presence of the exponent 7 = 1/3 seems to call for a different stochastic 
model. 

Finally, let us recall that the exponent 7 can provide information about the localization 
of Syl^ only when the convergence of the SLS is sufficiently fast. For instance, the exponent 
7 = 1/2 is observed also for a homogeneous chain, where no localization is obviously present. 
The apparent contradiction is solved by noticing that in this case the T-th spatial exponent 
is exactly equal to for any T. Therefore, 7 does not describe the convergence of the 
minimum positive exponent as expected, and the above arguments are no longer valid. 



D. Singularities in the TLS 

For the sake of completeness, we briefly discuss also an example of a singularity occurring 
in the TLS, which is not related to localization features of Lyapunov vectors, but rather to 
the existence of a conserved parameter. 

All singularities in SLS discussed in the previous section are characterized by 7 < 1, 
that is by a vanishing density of Lyapunov exponents. The study of TLS in models with 
a conserved quantity, such as Eq. (H), has revealed the existence of the complementary 
phenomenon too, i.e. a power-law divergence in the density of exponents 

dnx 1 



around A = |3^. By calling ric the integrated density determined by the condition 
A(0, ric) = 0, the above equation can be rephrased as |A| ~ \n\ — nc\" with u = 1/(1 — 13). 
Fig. 10 shows the scaling behaviour around the critical integrated density for model @ 
with g{x) = El sin(27ra;) + e2X {ei^2 are two coupling constants). Crosses and diamonds cor- 
respond to approaching A = from above and below, respectively. The different slopes 
about 1.5 in the former case and 2 in the latter) confirm the existence of two distinct critical 



behaviours as conjectured in |30|. This phenomenon appears to be very general, although a 



global explanation is still lacking. 
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VI. CONCLUDING REMARKS 



In this first part, we have introduced and utihzed a general representation for describing 
the hnear evolution of perturbations in spatially extended systems. We claim that spatial 
and temporal spectra allow for a comprehensive description of all phenomena occurring in 
this framework. In the second part wc shall show how the comoving and rotated exponents 
mentioned in the Introduction can be derived from SLS and TLS. At the present moment it is 
still unclear to what extent spatial and temporal exponents are independent of one another. 
As far as the border dT> in the (/x, A)-plane is concerned, the two classes are practically 
equivalent. However, no relation has yet been found between the corresponding densities. 
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FIGURES 



FIG. 1: Plot of the boundary dV for the logistic CML for two values of the coupling (a) 
e = 1/3, (b) e = 2/3. 

FIG. 2: (a-c) Temporal Lyapunov spectra for the logistic CML (e = 1/3) for /i = , 
1.31 (i.e. /ic) and 3.0, respectively; (d-f) Spatial Lyapunov spectra for the same diffusive 
coupling and A = —2.0 , and 2.0, respectively. 

FIG. 3: Schematic plot of the contour lines of A (solid line) and /i (dashed line) and the 
boundary dB in the (n^, nx)-p\ane . 

FIG. 4: Plot of the boundary &D for the asymmetric Bernoulli chain for e = 1/3 : (a) 
a = 3/4, (b) a = l. 

FIG. 5: Plot of the boundary dT> for a spatiotemporal periodic orbit (period 3 in space 
and 2 in time). 

FIG. 6: Plot of the boundary &D for a frozen random pattern. Dots and lines refer to 
temporal and spatial exponents respectively. Below /i = /i^, dVx 7^ dV^j^. 

FIG. 7: Plot of the boundary dV for coupled standard maps (a) k = 0.5 (b) k = 4.0. 
Symbols refer to spatial exponents while the line is obtained from the temporal spectrum. 

FIG. 8: Scaling of the maximum temporal exponent around /x = 0, logistic maps with 
random (circles) or deterministic (diamonds) dynamics, e = 1/3. The dashed line has slope 
2. 

FIG. 9: Singularity of the spatial spectra {e = 1/3): (a) skewed piecewise CML with 
deterministic (diamonds) and random (crosses) multipliers, cubic CML (squares) (1/2 sin- 
gularity); (b) logistic CML with deterministic (triangles) and random (plus) multipliers, 
symplectic CML (circles) (1/3 singularity). Details on the mapping functions and the best- 
fit values of the exponents are reported in Table 1. 

FIG. 10: Singularity of the temporal spectra of CML (^ with conserved quantity: g{x) = 
£1 sin(27ra;) + e2X, ei = 1/5 and £2 = 1/20. The critical value of the integrated density is 
= 0.643. 
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TABLES 



CML model 7 



tent 0.37 

tent (rnd) 0.36 

logistic 0.32 

logistic (rnd) 0.32 

linearly coupled 0.32 

symplectic 0.31 



skewed Bernoulli 0.54 
skewed Bernoulli (rnd) 0.45 
cubic 0.49 



TABLE I. Singularity exponent 7 of the SLS for various CML models Eq. (^) with the mapping 
function of the interval [0, 1] being f{x) = 1 — |2x — 1| (tent), f{x) = 4x(l — x) (logistic), /(x) = 3x 



for X < 1/3 and /(x) = (3x — l)/2 otherwise (skewed Bernoulli), /(x) = 3/2x + x — x^ (mod 1) 
(cubic). Linearly coupled maps refer to Eq. (0) with F{x) = 4x/(x^ + 1) and symplectic to Eq. 
with k = 4.0. In every case a coupling e = 1/3 is taken. The label "rnd" refers to random 
multipliers, with the same distribution as that of the corresponding CML. 
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